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Abstract 


Climate  variability  is  simply  represented  by  teleconnection  patterns  such  as  the 
Arctic  Oscillation  (AO),  Antarctic  Oscillation  (AAO),  North  Atlantic  Oscillation  (NAO), 
Pacific/North  American  Pattern  (PNA),  and  Southern  Oscillation  (SO)  with  associated 
indices.  Two  approaches  can  be  used  to  predict  the  indices:  forward  and  backward 
methods.  The  forward  method  is  commonly  used  to  predict  the  index  fluctuation  p  at 
time  t  with  a  given  temporal  increment  z  .  Using  this  method,  it  was  found  that  the  index 
(such  as  for  NAO)  has  the  Brownian  fluctuations.  On  the  base  of  the  first  passage  time 
(FPT)  concept,  the  backward  method  is  introduced  in  this  study  to  predict  the  typical 
time  span  (r)  needed  to  generate  a  fluctuation  in  the  index  of  a  given  increment  p  .  After 
the  five  monthly  indices  (AO,  AAO,  NAO,  PNA,  SO)  running  through  the  past  history, 
the  FPT  density  functions  are  obtained.  FPT  presents  a  new  way  to  detect  the  temporal 
variability  of  the  climate  indices.  The  basic  features  for  the  index  prediction  are  also 
discussed. 
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1.  Introduction 


Complexity  in  climate  systems  makes  the  prediction  difficult.  One  way  to 
simplify  the  climate  systems  is  to  represent  low-frequency  variability  of  atmospheric 
circulations  by  teleconnection  patterns,  such  as  the  Arctic  Oscillation  (AO),  Antarctic 
Oscillation  (AAO),  North  Atlantic  Oscillation  (NAO),  Pacific/North  American  Pattern 
(PNA),  and  Southern  Oscillation  (SO).  Temporally  varying  indices,  s(t),  were  calculated 
for  these  patterns.  Here,  t  denotes  time. 

Among  them,  the  SO  index  (SOI)  was  first  to  use  showing  equivalent  barotropic 
seasaw  in  atmospheric  pressure  between  the  southeastern  tropical  Pacific  and  the 
Australian-Indonesian  regions  (Walker  and  Bliss  1937).  A  popular  formula  for 
calculating  the  monthly  SOI  is  proposed  by  the  Australian  Bureau  of  Meteorology, 

s(t)  =  l0xPm(t)z{Pm)  (1) 

SD(pm) 

Here,  pdm  is  the  mean  sea  level  pressure  of  Tahiti  minus  that  of  Dawin  for  that  month; 
{Pdm )  's  the  long  term  average  of  pditr  for  the  month  in  question;  and  SD(/?diff)  is  the  long 

tenn  standard  deviation  of  pdm  for  the  month  in  question.  The  SOI  ranges  from  about  -35 
to  about  +35. 

The  indices  of  the  other  oscillations  (AO,  AAO,  NAO,  PNA)  are  the  expansions 
of  similar  seasaw  phenomenon  such  as  the  Northern  Hemisphere  annular  mode  [or  called 
AO]  and  the  two  patterns  (i.e.,  NAO  and  AO)  as  “two-paradigms  of  the  same 
phenomenon”  (Wallace  2000).  These  oscillation  patterns  are  usually  calculated  on  the 
base  of  the  Rotated  Principal  Component  Analysis  (RPCA)  (Barnston  and  Livezey  1987). 
Monthly  (AO,  AAO,  NAO,  PNA)  indices  are  constructed  by  projecting  the  monthly 
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mean  (1000-,  700-,  500-,  500-hPa)  height  anomalies  onto  the  leading  EOF  mode.  The 
time  series  are  normalized  by  the  standard  deviation  of  the  monthly  index.  Since  the 
loading  pattern  of  (AO,  AAO,  NAO,  PNA)  is  obtained  using  the  monthly  mean  height 
anomaly  dataset,  the  index  corresponding  to  each  loading  pattern  becomes  one  when  it  is 
nonnalized  by  the  standard  deviation  of  the  monthly  index.  Detailed  information  is 
referred  to  the  website:  http://www.cpc.noaa.gov/. 

Long-term  time  series  of  the  climate  indices  represent  invaluable  sources  of 
information  to  represent  the  complex  seasonal,  interannual,  and  interdecadal  variabilities 
of  the  climate  systems.  For  example,  sustained  negative  values  of  the  SOI  often  indicate 
El  Nino  episodes.  These  negative  values  are  usually  accompanied  by  sustained  wanning 
of  the  central  and  eastern  tropical  Pacific  Ocean,  a  decrease  in  the  strength  of  the  Pacific 
trade  winds,  and  a  reduction  in  rainfall  over  eastern  and  northern  Australia.  Positive 
values  of  the  SOI  are  associated  with  stronger  Pacific  trade  wind  and  warmer  sea 
temperatures  to  the  north  of  Australia,  popularly  known  as  a  La  Nina  episode.  Waters  in 
the  central  and  eastern  tropical  Pacific  Ocean  become  cooler  during  this  time.  Together 
these  give  an  increased  probability  that  eastern  and  northern  Australia  will  be  wetter  than 
normal.  The  SOI  has  been  used  in  numerous  studies  as  an  indicator  of  the  status  of  the 
ENSO  phenomenon  and  as  a  predictor  for  global  and  regional  climatic  prediction  such  as 
temperature  (e.g.,  Smith  and  Sterns,  1993),  precipitation  (e.g.,  Ropelewski  and  Halpert, 
1987),  and  agricultural  products  (e.g.,  Rimmington  and  Nicholls,  1993).  The  other 
indices  have  similar  features. 

How  to  predict  effectively  the  climate  indices  has  practical  significance  due  to 
their  connections  to  the  large  scale  atmospheric  circulations.  Usually,  these  indices  are 
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treated  as  time  series  and  statistical  predictions  are  conducted  (forward  method).  For 
example,  the  singular  spectrum  analysis  (Keppenne  and  Ghil  1992),  the  wavelet  analysis 
(Torrence  and  Campo  1998),  and  the  nonlinear  analog  analysis  (Drosdowsky  1994)  were 
used  to  obtain  the  dominant  frequencies  of  the  SOI  time  series.  The  power-law 
correlations  were  found  for  the  self-affine  properties  of  SOI  (Ausloos  and  Ivanova  2001). 

It  may  possible  use  an  alternative  “backward”  method,  which  predicts  a  typical 
time  span  ( r  )  needed  to  generate  a  fluctuation  in  the  index  of  a  given  increment  p  .  This 
method  is  established  on  the  base  of  the  first  passage  time  (FPT)  concept,  which  is  widely 
used  in  many  disciplines  such  as  physics,  chemistry,  biology,  and  economics,  etc.,  but 
not  in  meteorology  until  recently  when  Chu  et  al.  (2002a,  b)  used  FPT  to  study  the  ocean- 
atmospheric  model  predictability.  The  major  purpose  of  this  paper  is  to  show  the 
usefulness  of  FTP  on  the  climate  index  prediction. 

The  rest  of  the  paper  is  outlined  as  follows.  Section  2  discusses  the  two 
approaches  to  predict  the  climate  indices.  Section  3  shows  the  FPT  concept,  applied  in 
this  paper.  Section  4  investigates  the  sensitivity  of  the  FPT  density  function  to  the  index 
reduction  p .  Section  5  shows  the  power-law  dependence  of  optimal  FPT  on  the  index 
reduction  density  p  and  classifies  the  variability  of  the  climate  indices  as  the  Brownian 
motion  using  the  FPT.  Finally,  Section  6  presents  the  conclusions. 

2.  Forward  and  Backward  Approaches  of  Prediction 
Monthly  varying  climate  indices  from  the  NOAA  Climate  Prediction  Center  show 
randomness  (Figure  1)  with  poor  predictability  because  their  phases  and  amplitudes  are 
rather  unpredictable  with  both  involving  many  (time  and  space)  scales  which  are  often 
intrinsic  to  chaotic  behavior;  see  reviews  for  the  NAO  index  (Palmer  2000;  Greatbatch 
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2000;  Wanner  et  al.  2001).  The  histograms  of  the  monthly  mean  indices  show  Gaussian- 
type  distributions  with  mean  values  near-zero  for  all  the  indices  (AO:  -0.136,  AAO: 
-0.01,  NAO:  -0.004,  PNA:  0.018,  SOI:  -0.077)  and  with  standard  deviations  near  1  for 
AO  (0.981),  AAO  (0.989),  NAO  (0.985),  and  near  10  for  SOI  (10.5)  (Figure  2). 

Both  forward  and  backward  approaches  exist  in  index  prediction.  The  forward 
approach  predicts  the  change  of  the  index  ( p )  at  time  t  with  a  given  temporal 
increment  r  from  analyzing  single  or  multiple  time  series.  Due  to  stochastic  nature,  the 
probability  density  function  (PDF),  p[p{t),  r] ,  should  be  first  constructed.  Collette  and 

Ausloos  (2004)  analyzed  the  NAO  monthly  index  (single  time  series)  from  1825  till  2002 
and  found  that  the  long  range  time  correlations  are  similar  to  Brownian  fluctuations.  The 
distribution  functions  of  the  NAO  monthly  index  fluctuations  have  a  form  close  to  a 
Gaussian,  for  all  time  lags.  This  indicates  the  lack  of  predictive  power  of  the  present 
NAO  monthly  index.  Lind  et  al.  (2005)  used  the  standard  Markov  analysis  to  get  the 
Chapman-Kohnogorov  equation  for  the  conditional  PDF  of  the  increments  p  of  the 
NAO  index  over  different  time  intervals  r  and  to  compute  the  diffusion  and  drift 
coefficients  ( D{1\  D{2))  from  the  first  two  moments  of  such  probability  distribution.  The 
random  variable  p(t)  is  found  to  satisfy  the  Langevin  equation 

^  =  D{1)  [p{t),  t ]  +  7(0^(2)[p(0,f]  (2) 

where  r/(t )  is  a  fluctuating  8  -correlated  force  with  Gaussian  statistics,  i.e., 

{ri{t)ri{t'))  =  28it-t'). 

Maharaj  and  Wheeler  (2005)  predict  the  daily  bi-variate  index  (multiple  time  series)  of 
the  Madden-Julian  oscillation  using  seasonally  varying  vector  autoregressive  (VAR) 
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models.  The  first-order  VAR  model  on  the  original  (non-differenced)  time  series  was 
found  to  be  the  most  satisfactory  for  forecasting  the  index  beyond  a  few  days.  Although 
this  model  shows  no  strong  skill  advantage  over  a  lagged  regression  technique,  it  has  the 
convenience  of  employing  only  a  single  set  of  equations  to  make  predictions  for  multiple 
forecast  horizons. 

The  backward  approach  predicts  the  typical  time  span  (r)  needed  to  generate  a 
fluctuation  in  the  index  of  a  given  increment/? .  To  do  so,  the  FPT  concept  is  used.  As 
the  index  date  runs  through  the  past  history,  the  accumulated  values  of  the  FPTs  form  the 
PDF  (called  FPT  density  function).  The  mode  of  the  PDF  shows  the  most  probable  time 
when  the  index  decreases  with  the  value  of  p  at  the  first  time. 

The  FPT  is  important  from  a  prediction  point  of  view  in  several  ways.  First, 
although  the  climate  indices  are  related  to  certain  circulation  patterns  such  as  the  negative 
SOI  values  generally  connecting  to  the  El  Nino  events,  one  still  does  not  know  the  exact 
time  of  the  El  Nino  onset.  Therefore,  the  best  one  can  do,  from  a  statistical  point  of  view, 
is  to  make  a  prediction  at  a  time  that  is  probabilistically  favorable.  This  optimal  time,  as 
we  will  see,  is  detennined  by  the  maximum  of  the  FPT  density  function,  i.e.  the  optimal 
FPT.  Second,  but  not  least,  the  FPT  density  function  will  by  itself  give  invaluable,  non¬ 
trivial  information  about  the  stochasticity  of  the  climate  indices.  Third,  FPT  effectively 
represents  the  ocean-atmospheric  model  predictability  (Chu  et  ah,  2002a,  b). 

3.  FPT 

In  a  series  papers,  Chu  et  al.  (2002a)  introduced  the  FPT  ( r  )  concept  into  the 
ocean-atmospheric  model  predictability  with  the  model  error  (i.e.,  ‘some  quantity’)  first 
exceeding  a  predetermined  tolerance  level  (i.e.,  ‘predetermined  criterion’),  found  the  FPT 
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density  function  satisfying  the  backward  Fokker-Planck  equation,  and  obtained  the 
analytical  solution  of  the  FPT  density  function  for  the  nonlinear  simplified  low-order 
Lorenz  atmospheric  model  (Nicolis  1992).  Both  linear  and  nonlinear  perspectives  of 
forecast  errors  are  analytically  investigated  (Chu  et  al.  2002b).  Furthermore,  the  FPT  was 
also  used  to  evaluate  the  full  physical  nowcast/forecast  ocean  prediction  system.  For 
example,  the  FPT  density  function  is  asymmetric  with  a  broader  and  long  tail  in  higher 
value  side,  which  indicates  long-term  predictability.  The  long-term  (extreme  long) 
predictability  is  not  an  “outlier”  and  shares  the  same  statistical  properties  as  the  short¬ 
term  predictions  (Chu  et  al.  2002c).  The  FPT  is  also  used  to  verily  the  model  ability  to 
predict  Lagrangian  drifter  trajectories  (Chu  et  al.  2004)  and  the  regional  ocean  model 
predictability  to  stochastic  perturbations  in  initial  conditions,  winds,  and  open  boundary 
conditions  (Chu  et  al.  2005). 

Using  the  same  concept,  the  FPT  is  used  to  explore  the  statistical  features  of  the 
time  series  of  climate  indices.  They  change  either  positively  or  negatively  at  a  given  time 
(Figure  1).  Then,  of  course,  it  is  interested  in  predicting  the  exact  change  at  a  point  in 
time.  However,  it  is  not  possible.  Therefore,  the  best  one  can  do,  from  a  statistical  point 
of  view,  is  to  predict  the  time  that  is  probabilistically  favorable  for  the  given  index 
change.  This  optimal  time,  as  we  will  see,  is  determined  by  the  maximum  of  the  first 
passage  time  distribution,  i.e.  the  optimal  FPT. 

Given  a  fixed  value  of  an  index  reduction  (p),  the  corresponding  time  span 
(positive)  is  estimated  for  which  the  index  reduction 

/*(*)  =  *(?  + At) -s(t),  (3) 

reaches  the  level  p  for  the  first  time, 
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rp(t)  =  inf  {Af  >  0  |  /JO  <  ~p)  , 


(4) 


which  is  called  the  FPT.  It  is  noted  that  similar  definition  can  be  used  for  the  index 
enhancement.  As  the  index  data  run  through  the  past  history  which  is  represented  by  (2), 
the  FPT  density  function  satisfies  the  backward  Fokker-Planck  equation  (Chu  et  al. 
2002a) 

^-rD<V,ol^-V^m(/vA=o.  <5) 

dt  L  }dp  2  H  ’  dp2 

The  cumulative  distribution  is  introduced  for  the  transition  time  being  larger  than  r p  ,  i.e. 

P(jp)=  f  p(r)dT.  (6) 

Jtp 

All  the  five  climate  indices  show  similar  Gaussian- type  distributions  (Figure  2). 
The  Brownian  fluctuations  identified  in  NAO  monthly  index  (Collette  and  Ausloos  2004) 
may  extend  to  the  other  indices.  For  the  Brownian  fluctuations  with  zero  drift,  the  FPT 
density  function  has  the  analytical  solution  for  finite  p  (Rangarajan  and  Ding  2000), 

p(r)  =  -j=-^JTex p  ,  a>0,  r>0.  (7) 

yj7T  T  \  T  J 


Here,  a  is  a  function  of  p .  With  p—>  0  as  a  limit  case  (i.e.,  no  index  reduction),  the 
FPT  tends  to  0  with  a  probability  of  100%.  Thus,  the  FPT  density  function  cannot  be 
represented  by  (7)  but  the  5  function, 

p(z)  — >  S(j),  as  p  — >  0  . 

The  maximum  value  of  p{ r)  (for  finite  p  )  indicates  the  most  favorable  FPT, 


(8) 
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With  the  known  FPT  density  function  (7),  the  k- th  moment  (k  =  1,  2,  ...)  of  FPT  for 
finite  pis  calculated  by 

oo 

Mk{p)  =  k^p(r)TkAdz,  k  =  l,...,°o .  (9) 

0 

The  mean  and  variance  of  FPT  can  be  calculated  from  the  first  two  moments 

(r„)  =  M,,  (10) 

(< Sr;)  =  M1-M (11) 

where  the  bracket  denotes  the  average  over  realizations. 

4.  Sensitivity  of  p(  rp )  with  Respect  to  p 

The  theoretical  distribution  (7)  solely  depends  on  the  parameter  a.  This  parameter 
can  be  determined  from  the  monthly  index  data.  First,  the  FPT  density  functions  are 
constructed  and  fitted  to  (7)  for  each  index  with  various  values  of  index  reduction p . 
For  example,  Figure  3  shows  p{rp)  with  p  =  1.9  for  (AO,  AAO,  NAO,  PNA)  indices 
and  with  p  =  25  for  SOI.  Note  that  almost  an  order  of  difference  in  p  -values  (1.9  versus 
25)  is  selected.  This  is  due  to  the  different  ranges  of  index  fluctuation:  (-4,  4)  for  (AO, 
AAO,  NAO,  PNA)  indices  and  (-40,  40)  for  SOI.  It  exhibits  a  rather  well  defined  and 
pronounced  maximum,  followed  by  an  extended  tail  for  very  long  FPTs  indicating  a  non¬ 
zero  and  important  probability  of  large  passage  times  (note  that  the  x  -axis  is 

logarithmic).  These  long  (toward  El  Nino  for  SOI)  FPTs  reflect  periods  where  the 
tropical  Pacific  is  in  a  strong  El  Nino  phase  and  needs  a  long  period  of  time  before 
finally  coming  to  an  even  stronger  El  Nino.  The  short  FPTs  on  the  other  hand  -  those 


10 


around  the  maximum  -  are  in  La  Nina  periods,  which  appear  to  be  the  most  common 
scenario. 

Second,  for  each  index  reduction  ( p ),  the  parameter  (a)  is  determined  from  the 
Tp  data  fitted  to  (7).  From  the  scatter  diagram  for  each  index  (Figure  4),  the  linear 
relationship  between  ( a ,  p ), 

a  =  al+a2p,  (12) 

is  found  using  the  least  square  method.  Flere,  Eq.(12)  represents  the  (NAO,  NPA,  SOI) 
indices  quite  well,  but  not  too  well  for  the  (AO,  AAO)  indices  since  oscillations  around 
the  linear  trend  are  found.  Despite  oscillations,  the  linear  correlation  between  a  and  p  is 

very  evident  with  positive  regression  coefficient  a2  (Table  1). 

To  better  understand  the  tail  of  this  distribution,  various  values  of  p  are 
considered.  If  this  level  is  small  enough,  it  is  likely  that  the  index  reduction  will  break 
through  the  level  after  the  first  month,  while  larger  FPTs  will  become  more  and  more 
unlikely.  However,  the  probability  for  a  large  FPT  value  will  not  be  zero;  if,  say,  a  small 
level  p  ,  then  a  period  of  strong  El  Nino  will  result  in  a  r  (t)  that  might  be  considerably 

larger  than  one  month  since  it  takes  time  to  recover  from  the  El  Nino  event.  For  instance, 
after  the  1992  El  Nino  event,  it  took  five  years  to  reach  a  new  El  Nino  event  in  1997. 

5.  Characteristics  of  the  Index  Variability  -  Brownian  Motion 

A  random  process  is  called  fractional  Brownian  motion  if  its  cumulative  FPT 
density  function  satisfies  the  power  law  (Ding  and  Yang  1995), 

(13) 
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withO  <  H  <  1 .  Here  H  is  the  Hurst  exponent.  For  H  =  1/2,  the  random  process  is  the 
ordinary  Brownian  motion.  The  cumulative  distribution,  P(r  ),  for  each  climate  index 

data  (Figure  5)  shows  a  power  law  feature  in  the  tail  of  the  distribution  scales.  For  a  very 
small  value  of  the  index  reduction,  p0  (0.01  for  the  AO,  AAO,  NAO,  and  PNA  indices 

and  0.1  for  SOI),  the  cumulative  FPT  density  function  calculated  from  the  index  data 
shows  that 

P(T0)~T-0a°  ,  (14) 

with  a0  — 1/2.  Since  the  Hurst  exponent  of  an  ordinary  Brownian  motion  is  II  =  1/2,  the 

empirically  observed  scaling  (see  Figure  5)  is  a  consequence  of  the  (at  least  close  to) 
Brownian  motion  behavior  of  the  climate  indices.  This  argument  of  an  unbiased 
Brownian  motion  is  also  strengthen  by  observing  that  (see  Figure  5) 

P(To  =  1)  1/2.  (15) 

This  indicates  that  the  climate  index  has  50%  chance  of  increase  and  decrease  at  each 
time  step  (one  month). 

Figure  5  also  shows  the  cumulative  distributions  P(r  )  for  different  values  of  p  , 
i.e.,  p  =  0.5,  1.0,  and  1.5  for  the  (AO,  AAO,  NAO,  PNA)  indices  and  p  =  10,  20,  and 
30  for  SOI.  From  this  figure  it  is  seen  that  the  tail  exponent,  a  ,  is  rather  insensitive  to 
the  value  of  p  .  In  particular  one  finds  that  a p~  1/2  over  a  broad  range  of  values  for  p  ,  a 

value  that  is  consistent  with  the  Brownian  motion  hypothesis. 

Moreover,  it  is  noted  that  as  the  level  p  is  increased  from  zero,  the  optimal  FPT 

(r*max))  moves  away  from  r0  =1  and  toward  larger  values.  How  does  the  optimal  FPT 
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r'liax'  depend  on  the  level  p  for  large/??  This  dependence,  as  measured  from  the 

empirical  distribution  (i.e.,  the  histogram  of  FPT),  is  shown  in  Figure  6.  Intuitively,  it  is 
clear  that  the  optimal  FPT  will  increase  as  p  increases. 

The  power  law  is  found  for  the  optimal  FPT  versus  p  , 

z{;m)~py  for  large/?,  (16) 

with  y~2.0  for  the  AO,  AAO,  and  SO  indices  and  y  ~~1 .0  for  the  NAO  and  PNA 
indices  (see  Figure  7),  which  consists  with  the  results  for  a  Brownian  motion  with  the 
theoretical  FPT  density  function  [see  (7)  and  (8)].  Substituting  (12)  into  (8)  gives 

r‘max)  =  a\ p1  +2 axa2p  +  af  .  (17) 

The  exponent  of  the  power-law  depends  on  the  regression  coefficients  (ara2)  (see  Table 
-1).  Large  p  usually  means  much  larger  than  1  for  the  SOI,  and  slightly  larger  than  1  for 
the  AO,  AAO,  NAO,  and  PNA  indices.  For  the  SOI,  the  regression  coefficients  are  given 
by  ax  =0.165,  a2  =0.125  ,  the  leading  tenn  of  (17)  is  alp2 ,  and  therefore,  the  exponent 
of  the  power-law  is  around  2.  For  the  AO  and  AAO  indices,  a2  is  larger  thanc^ .  The 
leading  term  of  (17)  is  a;p 2  for  large  p ,  and  therefore,  the  exponent  of  the  power-law  is 
around  2.  However,  for  the  NAO  and  PNA  indices,  a,  is  smaller  than  al .  The  leading 
tenn  of  (17)  is  laxa2p  for  large/?,  and  therefore,  the  exponent  of  the  power-law  is 
around  1.  Not  surprisingly,  there  are  some  deviations  from  standard  theories  for  the 
variation  of  the  optimal  FPT  with  the  level  p  (i.e.,  the  straight  line).  Furthermore,  we 
have  checked,  and  found,  that  there  is  an  approximate  symmetry  under 

p^-p 
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for  the  FPT  density  function.  One  therefore  does  not  have  to  consider  index  enhancement 
explicitly. 

6.  Conclusions 

(1)  FPT  presents  a  new  way  to  detect  the  temporal  variability  of  the  climate 
indices.  It  predicts  a  typical  time  span  ( r )  needed  to  generate  an  index  reduction  of  a 
given  increment  p .  FPT  is  a  random  variable  and  whose  density  function  satisfies  the 
backward  Fokker-Planck  equation.  Solving  this  equation,  it  is  easy  to  obtain  the  ensemble 
mean  and  variance  of  the  FPT  of  the  climate  indices. 

(2)  FPTs  for  the  five  climate  indices  show  the  Brownian  fluctuations.  This  means 
that  the  early  results  on  the  Brownian  fluctuations  for  the  NAO  index  (e.g.,  Collette  and 
Ausloos  2004)  are  also  valid  for  the  other  indices  (AO,  AAO,  PNA,  and  SO).  With 
8  — »  0  as  a  limit  case  (i.e.,  no  index  reduction),  the  FPT  density  function  tends  to  the 
8  function. 

(3)  For  a  very  small  value  of  the  index  reduction,  p0  (0.01  for  the  AO,  AAO, 
NAO,  and  PNA  indices  and  0.1  for  SOI),  the  cumulative  FPT  density  function  shows  the 
power-law  dependence  on  t  with  the  exponent  approximately  -1/2.  Another  well  known 

method  to  check  power-law  dependence  is  examination  of  the  autocorrelation  function. 
This  also  confirms  that  the  climate  indices  have  the  Brownian-type  fluctuations. 

(4)  The  optimal  FPT  has  a  power-law  dependence  on  the  index  reduction  (p)  for 
large  p  with  the  exponent  (~2)  for  the  AO,  AAO  and  SO  indices  and  (~1)  for  the  NAO 
and  PNA  indices. 

(5)  The  FPT  density  functions  as  well  as  the  variation  of  the  optimal  FPT  can  be 
applied  if  one  wants  to  estimate  the  most  probable  time  period  needed  for  the  low- 
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frequency  atmospheric  circulation  pattern  to  sustain  if  a  prediction  aims  at  a  specific 
optimal  transition. 

(6)  The  FPT  analysis  can  be  applied  to  other  air-ocean  time  series. 
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Table  Captions 


Table  1.  Coefficients  (al,a2)  in  Eq.(12)  for  the  five  climate  indices. 
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Figure  Captions 


Figure  1.  The  monthly  mean  climate  index  data  obtained  from  NOAA  Climate 
Prediction  Center  for:  (a)  Arctic  Oscillation  (AO),  (b)  Antarctic  Oscillation  (AAO),  (c) 
North  Atlantic  Oscillation  (NAO),  (d)  Pacific/North  American  Pattern  (PNA),  and  (e) 
Southern  Oscillation  (SO).  Detailed  information  is  referred  to  the  website: 
http://www.cpc.noaa.gov/. 

Figure  2.  Histograms  of  the  climate  indices:  (a)  Arctic  Oscillation  (AO),  (b)  Antarctic 
Oscillation  (AAO),  (c)  North  Atlantic  Oscillation  (NAO),  (d)  Pacific/North  American 
Pattern  (PNA),  and  (e)  Southern  Oscillation  (SO). 

Figure  3.  FPT  density  functions  with  particular  index  reduction  p  for  (a)  Arctic 
Oscillation  (AO),  (b)  Antarctic  Oscillation  (AAO),  (c)  North  Atlantic  Oscillation  (NAO), 
(d)  Pacific/North  American  Pattern  (PNA),  and  (e)  Southern  Oscillation  (SO). 

Figure  4.  Linear  relation  ship  between  the  parameter  a  in  the  FPT  density  function  and 
the  index  reduction  8  for  (a)  Arctic  Oscillation  (AO),  (b)  Antarctic  Oscillation  (AAO), 
(c)  North  Atlantic  Oscillation  (NAO),  (d)  Pacific/North  American  Pattern  (PNA),  and  (e) 
Southern  Oscillation  (SO). 

Figure  5.  The  empirical  cumulative  density  functions,  P(rp ),  for  different  values  of  the 
index  reduction  p  for  (a)  Arctic  Oscillation  (AO),  (b)  Antarctic  Oscillation  (AAO),  (c) 
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North  Atlantic  Oscillation  (NAO),  (d)  Pacific/North  American  Pattern  (PNA),  and  (e) 
Southern  Oscillation  (SO).  It  shows  the  power-law  features. 

Figure  6.  FPT  density  function  of  SOI  with  different  values  of  p  :  (a)  15,  (b)  20,  (c)  25, 
and  (d)  30.  It  is  noted  that  the  optimal  FPT  increases  as  p  increases. 

Figure  7.  Dependence  of  the  optimal  FPT,  r"i:,x\  on  the  index  reduction  p  for  (a)  Arctic 

Oscillation  (AO),  (b)  Antarctic  Oscillation  (AAO),  (c)  North  Atlantic  Oscillation  (NAO), 
(d)  Pacific/North  American  Pattern  (PNA),  and  (e)  Southern  Oscillation  (SO).  Power-law 
exists  for  large  p  values. 
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Table  1.  Coefficients  (alta2)  in  Eq.(12)  for  the  five  climate  indices. 


AO 

AAO 

NAO 

PNA 

SO 

al 

0.430 

0.508 

0.604 

0.600 

0.165 

a2 

0.696 

0.562 

0.351 

0.373 

0.125 
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Figure  1.  The  monthly  mean  climate  index  data  obtained  from  NOAA  Climate 
Prediction  Center  for:  (a)  Arctic  Oscillation  (AO),  (b)  Antarctic  Oscillation  (AAO),  (c) 
North  Atlantic  Oscillation  (NAO),  (d)  Pacific/North  American  Pattern  (PNA),  and  (e) 
Southern  Oscillation  (SO).  Detailed  information  is  referred  to  the  website: 
http://www.cpc.noaa.gov/. 


24 


nao 


10 


2*5 


c 

(V 


cr 

d> 


(e) 


pna 


mean:  -0.0767 
STD:  10.5 


-nnnlT 


-60  -40 


-20 


20 


40 


soi 


Figure  2.  Histograms  of  the  climate  indices:  (a)  Arctic  Oscillation  (AO),  (b)  Antarctic 
Oscillation  (AAO),  (c)  North  Atlantic  Oscillation  (NAO),  (d)  Pacific/North  American 
Pattern  (PNA),  and  (e)  Southern  Oscillation  (SO). 
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(a)ao:p=-1.9  a=1.81  (b)  aao:  p=-1 .9  a=1.88 


x  [month]  x  [month] 

(e)  soi:  p=-25  a=1.18 
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Figure  3.  FPT  density  functions  with  particular  index  reduction  p  for  (a)  Arctic 
Oscillation  (AO),  (b)  Antarctic  Oscillation  (AAO),  (c)  North  Atlantic  Oscillation  (NAO), 
(d)  Pacific/North  American  Pattern  (PNA),  and  (e)  Southern  Oscillation  (SO). 
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(a)  ao:  a  =  0.43  +  0.695*p  (b)  aao:  a  =  0.508  +  0.562*p 


(c)  nao:  a  =  0.604  +  0.351  *p 


(d)  npa:  a  =  0.6  +  0.373*p 


(e)  soi:  a  =  0.165  +  0.125*p 


Figure  4.  Linear  relation  ship  between  the  parameter  a  in  the  FPT  density  function  and 
the  index  reduction  8  for  (a)  Arctic  Oscillation  (AO),  (b)  Antarctic  Oscillation  (AAO), 
(c)  North  Atlantic  Oscillation  (NAO),  (d)  Pacific/North  American  Pattern  (PNA),  and  (e) 
Southern  Oscillation  (SO). 
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(a):  ao 


(b):  aao 


Tp[month] 


Figure  5.  The  empirical  cumulative  density  functions,  P(t  ),  for  different  values  of  the 
index  reduction  p  for  (a)  Arctic  Oscillation  (AO),  (b)  Antarctic  Oscillation  (AAO),  (c) 
North  Atlantic  Oscillation  (NAO),  (d)  Pacific/North  American  Pattern  (PNA),  and  (e) 
Southern  Oscillation  (SO).  It  shows  the  power-law  features. 
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(e):  soi 


Figure  7.  Dependence  of  the  optimal  FPT,  z^max),  on  the  index  reduction  p  for  (a)  Arctic 

Oscillation  (AO),  (b)  Antarctic  Oscillation  (AAO),  (c)  North  Atlantic  Oscillation  (NAO), 
(d)  Pacific/North  American  Pattern  (PNA),  and  (e)  Southern  Oscillation  (SO).  Power-law 
exists  for  large  p  values. 


